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ABSTRACT 



Aims. 

Methods. We used the dual-band receiver GREAT on board the SOFIA airborne telescope to perform observations of the [C n] 158 jim 
fine-structure line at the postitions of two giant molecular clouds (GMC) in the center of IC 342 (GMCs C and E) and compared 
the spectra with corresponding ground-based data for low- and mid-/ CO and [C i]. We performed model calculations assuming a 
clumpy photo-dissociation region (PDR) environment using the KOSMA-t PDR model code to derive physical parameters of the 
local medium. 

Results. The [C n] 158 fim emission resembles the spectral signature of ground-based atomic and molecular lines, which indicates a 
common origin. The emission from GMC E can be decomposed into a cool, molecular component with weak far-ultraviolet (FUV) 
fields and low, mean densities of 10 3 cirr 3 and a strongly excited starburst/PDR region with higher densities of 10 4 cm" 3 and FUV 
intensities of 250-300 Draine fields. The emission from GMC C is consistent with gas densities of 5000 cirr 3 , FUV intensities of a 
few Draine fields and total gas masses of 20 x 10 6 M Q . 

Conclusions. The high spectral resolution of the GREAT receiver allowed us to decompose the [C n] emission of the GMC E into 
a strongly excited gas component resembling a PDR/starburst environment and a quieter, less excited gas component and to analyze 
the different components within a single beam individually. 

Key words, galactic: ISM — galactic: individual: IC 342 — radio lines: extragalactic — radio lines: ISM — atoms: [C II] 
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1. Introduction 

IC 342 is a gas-rich spiral galaxy with active star formation in 
its nucleus. IC 342 is located behind the Galactic plane und is 
therefore highly obscur ed; the distance to the nearly fa ce-on 
galaxy is still debated. Tikhonov & Galazutdinoval ((2010) give 
a distance of 3.9 + 0.1 Mpc derive d from stellar photom etry. 
Observatio ns of planetary nebulae dHerrmann et al . 2008) and 
cepheids dSaha et al.ll2002l) give distances of 3.4+0.2 Mpc. 

Within its central 30" two molecular arms of a mini-spiral 
end in a clumpy c entral ring of dense g as, which surrounds a 
young star cluster. iDownes et al. I (fl99l showed the presence 
of five giant molecular clouds, A to E, around the nucleus 
of IC 342 with masses of ~ 10 6 M . This structure can be 
seen in Fig. [1] which shows a map of line-integrated 12 CO(l- 
0) emission observed with BIMA (Berkeley Illinois Maryland 
Association). The central molecular ring surrounds a nuclear 
star cluster with active star formation and strong far-ultraviolett 
(FUV) radiation, illuminating the molecular ring and producing 
photo-dissociation regions (PDRs) on the side facing the central 
cluster. Energetiv FUV photons dissociate and ionize molecules 
and atoms in the gas and effectively heat the gas and dust via 
photoelectric heating. Consequently, PDRs strongly emit radia- 
tion from species that are abundant and ex cited under these con- 
ditions, like [C nl. lMeier & Tu rner (2005) find that the chemistry 



in the ring is a mixture of PDR gas and regions of denser, more 
shielded material. Especially the GMC complex B appears to 
show propertie s of PDRs produced by the surrounding starburst 
environment (IDownes et al.fl9 92: Israe l & Baasl2 003). The spa- 
tial sizes of the central GMCs and the infrared luminosity of 
the inner 400 pc of IC 342 are both similar to the center of our 
Galaxy. Especially the relative contribution of the diffuse mate- 
rial to the overall [C n] emission remains unknown but it appears 
that the PDRs remain largely confined to the central ring. 

Detailed knowledge of the H2 column densities, excitation 
temperatures, and densities in IC 342 exists from CO and its 
isoto pomers as a basis for comparison with t he \C u] lines 
(e.g. Ishizuki et al. 1990; Downes et al. 1992; Turner & Hurt 



1 19921: iWright et alj|1993|; iTurner et al.l ll993: 
Meier & Turner] |200irisrael & Baas 12003 



Meier et al. 2000; 
Meier & Turner! 



120051) . Comparison with this complementary data, particularly 
with their spectral shape, allows us to distinguish different kine- 
matic components from the [C 11] lines and to assess the interac- 
tion of the star formation activity on the gas. 

In this work we present velocity-resolved spectra of ion- 
ized atomic ([C 11] 158 fim) gas, with spatial resolutions of 
~15".6. The observations were performed with the German 
Receiver for Astronomy at Terahertz Frequencies^ (GREAT, 
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1 GREAT is a development by the MPI fiir Radioastronomie and 
KOSMA/ Universitat zu Koln, in cooperation with the MPI fur 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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Fig. 1. Line-integrated map of the 12 CO(1-0) transition from the 
BIMA-SONG sample. The black dots denote the postitions of our two 
[C n] detections in IC 342. The positions correspond to giant molecular 
clouds C and E. The (0,0) position corresponds to (RA,DEC) (J2000) 
(03:46:48.5 68:05:47). 



iHevminck et al.ll2012h on board the Stratospheric Observatory 
For Infrared Astronomy (SOFIA). 

In Sect. 2 we describe the observations. The spectra are pre- 
sented in Sect. 3. The analysis of the data and ambient conditions 
are discussed in Sect. 4, and conclusions are drawn in Sect. 5. 

2. Observations 

We used the dual-channel receiver GREAT on SOFIA in 
September 2011 to perform pointed observations of the 
[C n] 158 fan fine-structure transition of C + at 1900.536900 
GHz and the 12 CO(11-10) transition at 1496.922909 GHz near 
the center of IC 342. The observations were performed in dual 
beam-switch mode (chop rate 1 Hz) toward selected positions, 
see Tab.Q]for details. The chopper offset was 100" with a chop- 
per angle of 1 10 degrees against north (counter-clockwise). The 
observations were taken during the transfer flight from the US to 
Europe. Owing to unknown technical errors on the flight we had 
to discard one of the observed positions. 

The center position is RA,DEC (J2000) 03:46:48.5 68:05:47. 
We observed 2 GMCs at the following offsets in arcsecond: 
GMC C (+3.6,+3.6) and GM C E (-5.5,-3.7) using the designa- 
tion by iDownes et al. | ( [19921) (see also Fig. ffl). The integration 
times were 9.3, and 3.7 minutes, respectively. [C n] emission 
was detected at both positions. We did not detect any 12 CO(l 1- 
10) emission. The rms of the baseline for position C and E was 
0.10 and 0.21 Kkms~' , respectively, with a channel width of 4 
km s . 

We used a fast Fourier transform spectrometer (FFTS), with 
8 1 92 channels providing 1 .5 GHz bandwidth and about 212 kHz 
of spectral res olution. Ca libration was performed with the stan- 
dard pipeline dGuan et al. 2012). Using the beam efficiency tj c » 
0.51 and the forward efficiency (77/) of 0.95 (IHevminck et alJ 
1201 21) . we converted all data to line brightness temperature scale, 
Tb - r\f X T*Ji]c- The reduction of these calibrated data, as 
well as the maps shown throughout the paper, were made with 
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Fig. 2. Averaged [C 11] spectra (dark gray, filled) toward the GMC C 
(top), and GMC E (bottom) positions. The 12 CO(2-l) spectra are di- 
vided by 2, the [CI] 3 fi - 3 Pq and 13 CO spectra are multiplied with 
2. 



the GILDAS0 package CLASS90. The pointing was established 
with the optical guide cameras to about 1" precision, and was 
found to have deviated by about 3-5" after 20 minutes. 

3. Results 

In Fig. [2] we show the [C 11] emission observed with GREAT 
(dark gray, filled) for both positions. On e ach spectrum we over- 
lay 12 CO(1-0) data from BIMA-SONCfl H eifer et al.ll2003l) . as 
well as 12 CO(2-l), 12 CO(3-2), 12 CO(4-3 ), 13 CO(2-l) 1 3 CO(3 - 
2), and [CI] 3 Pi - 3 P spectra from llsrael & Baas! (120031) . 
Because of the different beam sizes, gridding, and map cover- 
age it was not always possible to regrid and convolve all spectra 
to the positions and beam sizes of the GREAT observations. 

Where possible, we smoothed the other data sets to the 15" 
resolution of our [C 11] line observations. We kept the 12 CO(2- 
1) and 13 CO(2-l) spectra on their native resolution of 22" and 
23". We kept 12 CO(3-2) and 13 CO(3-2) on their native resolu- 
tion of 15". 12 CO(1-0), 12 CO(4-3), and [C i] 3 Pi- 3 P spectra were 
smoothed to the resolution of our [C 11] data and interpolated at 
the two positions C and E. For the remaining ground-based spec- 
tra at the positions C and E we averaged the data points closest 
to our [C 11] observations. 

The [C 11] emission is strongest at position E. At both po- 
sitions the different transitions show a good correlation of the 
central line velocities and line widths. The spectra at position C 
exhibit a consistent central velocity of 40 kms~' and similiar 
line widths of about 50 km s . Only the [C 11] line has a broader 
line width and is blueshifted by a few km s -1 . Position E shows 
slight variations in the different spectral lines. All lines have a 
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comparable line width of slightly less than 50 km s . Most lines 
appear to be centered around 24 kms -1 , except for 12 CO(2-l) 
and (3-2) with slightly redshifted central velocities. 

We performed a Gaussian fit for all lines at the two posi- 
tions to determine their respective line parameters such as peak- 
and line-integrated intensity T^ k and J T m bdv, line width A V 
, and line center velocity vo. For position C we fitted a single- 
component Gaussian to the data without any additional con- 
straints. The results are given in Table Q] At position E we al- 
ready noted the slight asymmetry in several of the CO lines, 
possibly resulting from two different, spatially and kinematically 
unrelated gas components. We us ed the optically th innest lines 
from our data set, the 13 CO lines dlsrael & Baasl2003l) . to specify 
the line center velocity of 24 kms~' of the main component by 
fitting a single Gaussian to their profiles. To quantify the contri- 
butions of two separate components to the spectra, we fitted two 
separate Gaussian components to the spectra while keeping the 
line center velocity of the main component fixed to 24 kms~'. 
Only for 12 CO(4-3) and [C i] it was not possible to find a sec- 
ond line-component, therefore we performed an unconstrained 
single-component fit that we can attribute to the main compo- 
nent because of the comparable velocities. For 12 CO(1-0) we 
fixed the velocity and width of the high-velocity component to 
53 and 47 kms -1 , respectively, and derived the line parame- 
ters for the low-velocity component. The decomposition into two 
different components is also supported by the BIMA 12 CO(l- 
0) data at 5.5" resolution, which shows a line center velocity 
vo ~ 50 km s at our center position. Therefore, emission from 
the PDR-dominated gas close to the central cluster is expected 
at vo ~ 50kms~'. 

As an independant method to derive the [C n] line parame- 
ters, we also calculated the zeroth, first and second moment of 
the spectra, i.e. the integrated intensity, mean velocity and line 
width. The derived emission line parameters are summarized in 
Table Q] (moments are given as italic numbers). 

The GMC E and B are both located southwest of the cen- 
tral position of IC 342 with an angular separation of only 3.5". 
Accordingly, the 15" [C n] beam will always pick up emis- 
sion of both GMCs. The [C n] spectrum of position E shows 
the strongest emission from both positions most likely because 
of the already known P DR/starburst contribution from GMC B 
(Meier & Turner 2005). Position C is situated in the molecular 
ring, closest to the central cluster and will most likely contain 
contributions from dense molecular gas as well as from PDRs. 

Ilsrael & Baasl d2003l) showed the existence of strong veloc- 
ity gradients in the central region of IC 342. To quantify what 
pointing shift could alternatively mimic the observed line shift, 
we convolved the BIMA 12 CO (1-0) data to 15". We found that 
a pointing offset of 7 ", significantly higher than the quoted ac- 
cumulative pointing accuracy after 20 minutes, exactly in the 
direction of the central cluster could explain the [C n] line pro- 
file. Therefore, it is unlikely, though it cannot be fully excluded, 
that a systematic pointing error may have caused the observed 
velocity shift. 



4. Discussion 

Table |2] lists line ratios of [C n] and [C i] to 12 CO(1-0) and 
12 CO(2-l) and 12 CO(4-3)/ 12 CO(l-0) (if available, the numbers 
for position E are given in the order high/low velocity compo- 
nent). Stronger PDR emission will be reflect in a higher ratio 
of the fine structure lines to the molecular lines. Apparently, 
GMC E shows the strongest PDR contribution. IStacev et al.l 



Table 1. IC 342 line parameters derived from Gaussian fits. One 
Gaussian component was fitted to all spectra at position C. At 
position E we fitted two Gaussian components to the spectra with 
a line center velocity of component 1 fixed to 24 kms~' as de- 
rived from the single-component fit to the 13 CO(2-l) spectrum. 
The line parameters for [C n] are additionally calculated from 
their moments between -30 and 100 kms _1 (numbers in italics). 



Transition 


/ T„ b dv 


vo 


AV 


L mb 




[Kkms- 1 ] 


[tor 1 ] 


[kms" 1 ] 


[K] 


Position C (+3.6", +3.6") 


[C U] 


24.3±2.7 


33.1+4.5 


78+9 


0.3 


[C II] « 


23.5 


31.8 


67.5 




12 CO(1-0) 6 


307±2 


41.5+0.2 


50+1 


5.83 


12 CO(2-l) 


131±3 


38.2+0.5 


55+1 


2.21 


12 CO(3-2) 


134±1 


39.4+0.2 


53+1 


2.37 


12 CO(4-3) 


98±3 


39.7+0.9 


54+2 


1.71 


[CI] 3 P,- 3 P„ 


25.6±2.2 


38.3+2.0 


47+5 


0.51 


Position E (-5.5", -3.7") 


[C II] 


37±15 


24 r 


48+9 


0.72 




20.6±9.9 


52.7+11.3 


47+20 


0.42 


[C II] » 


59 


34.9 


51 




"coa-of 


280±11 


20 


43+2 


6.29 




21.8+9.4 


53 d 


47 f 


1.45 


12 CO(2-l) 


139+6 


24 c 


50+2 


2.59 




9.5+4.2 


59.8+4.5 


34+11 


0.26 


12 CO(3-2) 


120+5 


24 r 


49+1 


2.27 




16+3 


53.3+4.3 


46+8 


0.32 


13 CO(2-l) 


23.6+2.7 


24.0+1.7 


48+4 


0.46 


13 CO(3-2) 


13.9+2.0 


24.4+1.8 


41+4 


0.32 


12 CO(4-3) 


88+5 


21.2+0.9 


46+2 


1.81 


[C I] 3 P,- 3 P 


34±4 


24.4+1.7 


51+4 


0.61 



Notes. (fl) Results from calculating the moments of the line. (6) To 
convert from flux density Jy beanr'km s -1 to brightness temperatures 
we applied a facto r of 4.4 assuming an effective beam size of 5.5" 
(Heifer et al. 2003). (c) The line center velocity was set to 24kms~'. 
(rf) The line center velocity was set to 53km s . (e) The line width was 
set to 47km s -1 . 

(119911) give a [C n]/ I2 CO(l-0) ratio of 5000 (corrected for main 
beam efficiency and assuming unity beam filling). They cite 
[C n] intensities (beamsize 55") of 3xl0~ 4 erg s _1 crrr 2 sr _1 , 
but only 40 Kkms~' for the 12 CO(1-0) line intensity in a 
60" beam. The [C n] and 12 CO(1-0) data were taken using the 
Kuiper Airborn Observatory (KAO) and the Owens Valley Radio 
Observatory (OVRO), respectively. Given the larger beam sizes 
and unknown filling factors, their values are consistent with our 

data^ 

IStacev et al.1 dl99ll) showed that starburst galaxies have a 
high [C n]/ 12 CO(l-0) ratio of 4100 while cooler, less ac- 
tive galaxies show much weaker ratios. Decomposition of the 
I2 CO(1-0) data into a low- and high-velocity component at 24 
and 53 km s , respectively, gives [C n]/ 12 CO(l-0) ratios of 600 
and 4000, indicating contributions from quiescent, weakly ex- 
cited gas plus a strong, starburst like contribution from the local 
PDRs. 

The [C n] emitting level lies 91 K above the ground state and 
has a critical density of > 3500 cm 4 . Assuming optically thin 
emission and a level population in the high-temperature, high- 
density limit, we can calculate a lo wer limit of the C + co lumn 
density from the integrated intensity (Crawf ord et al.M l985). The 
results are given in Table [2] The beam-averaged column densi- 
ties have values between 0.9 - 1.6 x 10 17 cm" 2 , relatively low 
numbers for massive PDRs. The source intrinsic column density 
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Table 2. Line ratios and estimated column densities of the ob- 
served GMCs. For GMC E we give the values for the high/low 
velocity components. The PDR fit assumes a distance of 3.9 
Mpc. 





GMC C 


GMC E 


[C n]/ 12 CO(l-0) 


482 


4236/692 


[C n]/ 12 CO(2-l) 


142 


1223/150 


[C i]/' 2 CO(l-0) 


6.5 


-/9.4 


[C i]/' 2 CO(2-l) 


1.9 


-/2.4 


12 CO(4-3)/ 12 CO(l-0) 


20.7 


21.1(-) 


%:„]/[ 10 17 cm- 2 f 


1.4 


0.9/1.6 


PDR model results 


(n) [10 3 cm" 3 ] 


5.0 


10/2.0 


M tot [10 6 M G ] 


20 


2.0/15 


X [Draine] 


7 


300/5 



Notes. The line ratios were calculated from line-integrated intensities 
in units of erg s cirT 2 sr -1 . 

(o) Assuming a high-temperature and high-density limit and optical thin 
emission. 

will be accordingly higher for a compact source, lower than the 
beam: a beam filling factor of 1/10, not unlikely given the small 
size of the cores visible in the interferometric 12 CO(1-0) map, 
will give a column density that is a factor of 10 higher. Column 
densities of 1 - 2 x 10 17 cm" 2 and a [C + ]/[H] ratio of 10 -4 gives 
a total ma ss in the 15" beam o f 5 - 10 x 10 6 M Q at a distance of 
3.9 Mpc. Ilsrael & Baasl d2003l) give comparable central masses 
of IC 342 of 5 - 7 x 1 6 M , given the uncertainties in the relative 
abundance of ionized carbon, filling factors, and distance. 

To confine the lo cal gas parameters we used the KOSM A- 
t PDR model code dStorzer et al.lll996t iRollig et al.1 120061) to 
model t he emission of an e nsemble of clumpy PDRs in a beam 
of 15" (ICubick et al.ll2008l) . We fitted the model to the absolute 
intensities therefore making the total model gas mass sensitive to 
the distance of the PDR. As a result from the fit we receive the 
mean gas density in the beam, the total gas mass, and the FUV 
intensity, that illuminates the PDRs in units of the Draine field 
(|Prainell978l) . The derived parameters for all positions are given 
in Tab. [2] assuming a distance of 3.9 Mpc. Employing a shorter 
distance will reduce the total gas mass estimates accordingly. 

To perform a consistent fit across both positions we excluded 
the 13 CO lines from the fitting. Note that the I3 CO lines show 
relative weak emission, which requires lower model densities of 
~ 1000 crrT 3 for both positions. Without fitting to the 13 CO lines 
the densities are higher by a factor of a few. Furthermore, the 
[C n] and [C i] line intensities relative to the CO emission can 
only be explained by densities below 10 4 cm 3 and low FUV 
fields. Not e that the model co mputations assume solar metallic - 
ities, while lEngelbracht et al.l d2008l) give somewhat higher val- 
ues for IC 342. Higher metallicities will lead to lower [C n]/CO 
ratios. 

The fit results from GMC E support the decomposition into 
two separate gas components. The lower velocity of 24 kms 1 
would be consistent with assuming that the emission originates 
in the trailing, spiral arm region, which is connected to the 
molecular ring. The emission of the high-velocity component 
resembles the PDR signature of a starburst/PDR environment 
with mean FUV fields of 250-300 (in units of the Draine field) 
and densities of 10 4 crrT 3 . If we assume that the observations 
at position E were mispointed (see discussion in the previous 
section), we derive the following PDR model parameters: (n) = 



1000 cirT 3 , M to t = 19 x 10 6 M and^ = 5. These gas parameters 
are in conflict with a scenario where a PDR/starburst dominates 
at the center of the galaxy. 

5. Conclusions 

We used the dual-band heterodyne receiver GREAT on board 
the airborne telescope SOFIA to observe two giant molecular 
clouds situated around the nucleus of IC 342 in the I2 CO J = 
11 — > 10 transition and the [C n] 158 fim fine-structure line. We 
detected [C n] emission at both positions but could not detect 
any 12 CO J — 1 1 — > 10 emission. 

The new SOFIA/GREAT spectra reveal a spectral distribu- 
tion of the [C n] emission that follows the distribution of the 
neutral and molecular gas, [C i] and CO. The [C n] spectrum 
observed at the position of GMC E shows two velocity compo- 
nents, a high-velocity component that we attribute to emission 
from a PDR/starburst region in the molecular ring close to the 
central cluster with densities of 10 4 cm -3 , FUV field of 250-300 
and a total mass of 2 x 10 6 M and a cooler, low-velocity com- 
ponent with densities of 2 x 10 3 cirT 3 , FUV fields of a few and 
a total mass 7-8 times higher than the starburst component. The 
model for GMC C gives model parameters of densities of 5 x 10 3 
cirT 3 , FUV field of 7 and a total mass of 2 x 10 7 M . 

Despite the challenges that one might expect in an early 
transfer flight of SOFIA to Germany, we were able to deduce 
important astrophysical results primarily owing to the high spec- 
tral resolution available with the GREAT receiver. These data 
demonstrate the promise of the GREAT/SOFIA facility for fu- 
ture work, such as detailed mapping of the [C n] emission from 
the central regions of IC 342 and other nearby galaxies. 
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